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We propose a degree-based coarse graining approach that not just accelerates the evaluation of 
dynamics on complex networks, but also satisfies the consistency conditions for both equilibrium 
statistical distributions and nonequilibrium dynamical fiows. For the Ising model and Susceptible- 
Infected-Susceptible epidemic model, we introduce these required conditions explicitly and further 
prove that they are satisfied by our coarse-grained network construction within the annealed network 
approximation. Finally, we numerically show that the phase transitions and fluctuations on the 
coarse-grained network are all in good agreements with those on the original one. 
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I. INTRODUCTION 

Complex network has been one of the most active re- 
search topics in statistical physics and many other disci- 
plines [il-Q- It describes not only the pattern discovered 
ubiquitously in real world but also a unified theoretical 
framework to understand the inherent complexity in na- 
ture. However, the investigation of large networks, such 
as human brain that composes of about 10^^ neurons and 
10^ synapses [^], requires tremendous time-demanding 
efforts. The phenomenological description may capture 
certain properties of system, but always neglects micro- 
scopic information. A promising alternative is to develop 
coarse-grained (CG) methods, aiming at significant re- 
ducing the degree of freedom while proper preserving the 
microscopic information of interest. 

Several CG schemes have been proposed. Renormal- 
ization transformation has been used to simplify self- 
similar networks, and the reduced networks often pre- 
serve some topological properties of the original ones [Tl- 
[ig |. Spectral coarse graining technique has been pro- 
posed, in which the eigenvalues of Laplace matrix of 
network are almost unchanged, such that the dynam- 
ics of random walk and synchronization are preserved 
[ill [l3 | . Equation-free multiscale computational frame- 
work has been developed to accelerate simulation using a 
coarse time stepper [13]. This approach has been applied 
to study the CG dynamics of oscillators network _14] . 
gene regulatory network [15| , and adaptive epidemic net- 
work Il6|. However, to our best knowledge, no attempt 
has been made for developing CG simulation method to 
study critical phenomena (usually described by stochas- 
tic models) on complex networks. The size-dependent 
and scaling behaviors in these systems are studied so far 
mainly by such as Monte Carlo (MC) and kinetic MC 
(KMC) simulations [13, [l3- Apparently, these micro- 
scopic approaches are often too expensive. It is noticed 
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that the CG stochastic models have been proposed to 
study reaction-diffusion processes on regular lattices ^M, ■ 
However, the existing methods are largely inapplicable to 
critical phenomena on complex networks with diversified 
heterogeneity. Moreover, the crucial issue concerning the 
methodology development as to what criteria should be 
met to make the CG model statistically consistent with 
the microscopic one is yet to be addressed. 

In this paper, we propose the degree-based CG (d-CG 
for short hereafter) be a statistically consistent scheme, 
within the annealed network approximation (ANA) [iTI . 
[20 . [2H . It may therefore be an efficient and reliable CG 
method for evaluating the stationary properties of phase 
transitions and studying size effect on complex networks. 
We put forward the conditions of statistical consistency 
(CSC) on both equilibrium and nonequilibrium proper- 
ties, exemplified with the Ising model and Susceptible- 
Infected-Susceptible (SIS) model, respectively. We show 
that the d-CG approach that merges together the nodes 
of similar degrees does warrant the CSC within ANA. 
The stochastic Ising spin-flip and epidemic spreading 
dynamics can therefore be evaluated faithfully and ef- 
ficiently with the d-CG networks. The calculated phase 
diagrams, fluctuation dynamics, and system size-scaling 
behaviors are all shown in excellent agreements with the 
corresponding microscopic MC and KMC results. 

The rest of this paper is organized as follows. In SeclHl 
we present a general scheme for coarse graining network 
and d-CG approach, and further prove that this ap- 
proach satisfies statistical consistency. Extensive numer- 
ical demonstrations of the CG approach are performed 
on diverse networks in Sec lIIII At last, main conclusions 
and discussion are addressed in Sec lIVl 



II. COARSE GRAINING PROCEDURE 
A. Netw^ork Coarse Graining 

Let us start with the basic ingredients of network 
coarse graining. Consider a network consisted of N nodes 
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FIG. 1: (color online) (a) A schematic example of coarse 
graining network in which nodes of the same shape are 
merged; (b) the adjacency matrix A"^ of CG network by 
Eq. IT]). Note that the CG network may include self- 
connections; for example, the CG-node III where the con- 
nected nodes are merged. 



whose connectivity is given by the adjacency matrix A, 



in which Ai 



1 if nodes i and j are connected and 



Aij = otherwise. Merging q^ nodes together into a 
CG-node (denoted by C^) leads to a CG-network with 
N'^ CG-nodes. We adopt the mean-field definition of 
the CG connectivity between C^ and C,y, as the average 
number of links connecting any two nodes inside C^ and 
Cj/. The adjacency matrix A'^ of the CG network is then 



-^—^j: A,, \iii = v. 
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For illustration, Fig[T]depicts an example of coarse grain- 
ing a network with six nodes into a weighted CG network 
with three CG-nodes. 



B. CSC 

We address the issue on the statistical consistency of 
a CG scheme with the microscopic network, in terms of 
both the equilibrium distribution and the nonequilibrium 
fiow. We exploit the Ising model and the SIS model as 
the paradigmatic examples for equilibrium and nonequi- 
librium systems, respectively. 

For the Ising model defined on a network, the Hamil- 
tonian is given hy H = —J'^^^j AijSiSj, with the spin 
variable Si = ±1, and the ferromagnetic interaction pa- 
rameter J > 0. The probability of a given microscopic 
spin configuration {sj} is given by canonical distribution 
e~^^ /Z, where /? is the inverse temperature and Z is 
the partition function. The corresponding CG Hamilto- 
nian assumes the sum of pairwise interactions inside and 



between the CG-nodes, i.e. H'^ — H^ -t- iff, 
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Here, 

^± - 



r]^ = EiGC ^i ^^ *^^ ^^ ^P^'^ variable, and 
{q^ ± ri^)/2 is the number of up/down spins in- 
side Cfj_. Note that H'^ = iJf + iJ^ above is a closure 
expression at the CG level that depends only on A'^ and 
{ry^}. To make the CG model consistent with the micro- 
scopic one, it demands that the probability of any given 
CG configuration {77^} in equilibrium be the sum of the 
probabilities of all microscopic configurations that con- 
tribute to it. That is 
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where 5({?7/j}) = Hu ^^■/("m '"m ') ^^ ^^® number of mi- 
croscopic configurations corresponding to {?7/i}. 

We now consider the nonequilibrium scenario, exempli- 
fied by the SIS network for epidemic spreading dynamics. 
At the microscopic level, the SIS network nodes represent 
individuals being either susceptible {a = 0) or infected 
((7 = 1). A susceptible individual j can get infected at 
rate fj = Y^ih^i^ where fi^j = \Aij(Ti{l - aj) is the 
spreading fiow from an infected individual i to j. On the 
other hand, an infected node can recover to susceptible 
state at rate r^ = 70"^. Without loss of generality, we set 
7 = 1 hereafter to scale the infection parameter A. The 
SIS model exhibits a nonequilibrium dynamical phase 
transition at A = Ac from an absorbing state (all recov- 
ered) to an active state (disease spreading persistently) 
[H [23]. For the CG-SIS network, we define the CG vari- 
ables as the number of infected individuals inside a CG- 
node C^ by af^ = '^i^c '^i- Written in a closure form, 



the CG recovery rate is 



a'l^ and the CG spreading 



flow between two CG-nodes is J^^^, = AA^^cr^((7y — al). 
For this nonequilibrium system, the CG model is consis- 
tent with the microscopic one if the CG flow matches the 
microscopic flows: 
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Note that the CG recovery rate of r'i 
trivially. 



(5) 



T.i<^c, ""-^ holds 



C. Degree-based CG Scheme 

In our approach, we merge the nodes with similar de- 
grees together. We shall show the resulting d-CG scheme 
does satisfy the CSC, as defined by Eq. ^ for the Ising 
model and Eq. ([5]) for the SIS model, within the ANA for 
the ensemble averaged dynamics [l3, [20, [21|. In many 
previous studies [23!, .Mj , ANA has been extensively con- 
firmed to be a useful tool for describing quenched net- 
works, as in case of the present paper. Although ANA 
is just an approximation, it still gives a reasonable de- 
scription of the average behavior of nodes of the same 
degree. Moreover, one can consider that the ANA is a 
statistical characterization of large number of quenched 
networks with the same degree distributions. Neverthe- 
less, in Ref.[25| it has been pointed out that there exists 
some discrepancies between considering annealed approx- 
imation for quenched networks and considering annealed 
network models by themselves. According to the ANA, 
one can replace the dynamics on a given network of N 
nodes by that on a weighted graph of the full connectiv- 
ity Aij — kikj/{N{k)), where ki and kj are the degrees of 
nodes i and j, respectively, and (fc) is the mean degree. In 
an ideal d-CG scheme, the microscopic nodes in a single 
CG-node are of the same degree: ki\i^c —K^. The CG 
connectivity is then AJ^^ = IMW^ ^*^^ ^^^ CG-network 
dynamics treated at the ANA level. As results, Eqs.([2]) 
and ([3]) become, respectively (noting that s^ — 1) 
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Their sum can be written as 
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It is identical to the microscopic Hamiltonian at the ANA 
level. In other words, the Hamiltonian of any d-CG con- 
figuration equals to the collection of its contributing mi- 
croscopic configurations. The pre-exponential terms in 
two sides of Eq.(j4]) are both the degeneracy of CG config- 
uration, as well as the identical energy factors {H'^ — H) . 
We have thus proved that the ideal d-CG approach to the 
Ising model obeys the CSC of Eq.Q exactly. It is also 
easy to prove that the nonequilibrium CSC of Eq.([5|) is 
true for the SIS model in consideration; both sides there 
equal to A j^/^x o-f^iq^ ~ o-f,). Therefore, the d-CG ap- 
proach satisfies the CSC for both equilibrium probability 
distributions and nonequilibrium dynamical flows. Cer- 
tainly it is anticipated that the CSC holds approximately 
if merged together are the nodes of similar degrees rather 
than the exactly same ones. 
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FIG. 2: Typical time evolutions of the magnetization m in 
Ising model at T = 10 (in unit of J/ks) and the density of 
infected nodes p in SIS model at A = 0.1 for both microscopic 
and CG levels. Other parameters are N — 1024 and A'''^ — 16. 



III. NUMERICAL DEMONSTRATIONS 

A. CG-MC and CG-KMC simulations 

The MC simulation with the Metropolis dynamics [23| 
and the KMC simulation [23] are applied to the Ising 
model and SIS model, respectively, at both the mi- 
croscopic and the CG levels. In the CG-MC simula- 
tion, a CG-node C^ is randomly chosen, followed by the 
Metropolis try for spin-flip process. The probabilities 
of flipping a up/down-spin is n^W{l3,AE-f/^), where 
W{l3,AE) = min(e~'^'^-^, 1), and AE-^^^ is the change 
of CG Hamiltonian resulting from the flip of a up/down 
spin. It is easy to conflrm that the CG-MC simulation 
obeys the detailed balance condition. In the CG-KMC 
simulation, following the Gillespie algorithm |,27], a CG 
process to be executed is randomly selected based on the 
transition rates of all processes. The conflguration and 
transition rates are updated for executing the next CG 
process. 



B. Scale- free netw^orks 

We first consider the Barabasi-Albert (BA) scale- 
free network [28[, with the degree distribution follows 
a power-law P{k) ~ k^^ with scaling exponent 7 = 3. 
Figure [2] plots typical time evolutions of the magnetiza- 
tion m = ^^ Si/N in Ising model at T = 10 (in unit of 
J/ks) and the density of infected nodes p = ^^ (Ti/N in 
SIS model at A = 0.1, where N = 1024 and N" = 16 are 
used. For both the microscopic and CG simulations, the 
systems attain the steady states associated with fluctuat- 
ing noise after transient time. It is clear that there are in 
good agreement in the steady-state values of m and p, as 
well as their fluctuating amplitudes for both simulations 
cases. 
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FIG. 3: (a)-(b) \m\ and x as a function of T. Symbols and 
lines correspond to results by simulating the CG model and 
by direct microscopic MG simulation, respectively, (c) Tc as 
a function of In A^ obtained by the GG model with N'^ = 64. 
All dotted lines are the results of randomly merging models 
with iV = 16384 and N" = 64. All the networks have fixed 
mean degree (k) — 20. The error bars (not shown) are smaller 
than the symbol size. 
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FIG. 4: (a)-(b) p and and Ap as a function of A. Symbols 
and lines correspond to results by simulating the GG model 
and by direct microscopic KMG simulation, respectively, (c) 
Ac as a function of 1/ In A'^ obtained by the GG model with 
A^'^ = 64. All dotted line are the results of randomly merging 
models with iV = 16384 and A'''' = 64. All networks have 
fixed mean degree (fc) = 20. The error bars (not shown) are 
smaller than the symbol size. 



For the Ising model, we construct initial configurations 
by preparing each node with a random spin value Si = +1 
or —1 with an equal probability. As the simulation pro- 
ceeds, the system quickly relaxes to an equilibrium state. 
With the same temperature T we run at least 100 times 
of simulations corresponding to different initial configura- 
tions and network realizations. In each simulation, 2 x lO'^ 
MC steps (MCS: each spin is attempted to flip once on 
average during each MCS) are performed and the last 10'^ 
MCS are used to investigate the system's behavior. As 
T decreases the value of the magnetization m undergoes 
a transition from zero to nonzero at the critical tempera- 
ture T^. Below Tc we notice that due to finite-size effects 
the system can switch between two stable states via a nu- 
cleation mechanism, resulting in the oscillations of m 29] . 
Above Tc^ m fluctuates around zero, and the susceptibil- 
ity X per node has a maximum at phase transition, which 
can be used to determine Tc as we shall show in FiglHJb) . 
The susceptibility is related to the magnetization fluctu- 
ation via the fluctuation-dissipation theorem. To avoid 
the offset of m due to its oscillations in simulations, we 
will use instead the absolute value of m in the following. 
mi and X ^-re shown in Figj31Ja) and (b), respectively, 
as a function of the temperature T for different iV and 



N'^. Apparently, the d-CG results (reported by symbols) 
are in excellent agreements with the microscopic-level 
counterparts (by solid and dashed lines). As contrast, 
we also report the results of random-merging (RM) CG 
scheme (by dotted lines) in Figl^Ja) and (b), where each 
CG~node includes g^ = N/N" {N = 1024, N" = 64) 
nodes selected randomly from the whole network. Evi- 
dently the random scheme fails badly in reproducing the 
microscopic behaviors. We also used the same q^ in the 
random scheme as in the case of d-CG scheme, and found 
that the two random schemes produce the same results. 

It is important to merge nodes with similar degrees 
together, as already shown both analytically and numer- 
ically. Strikingly, even when the original network is re- 
duced to one with only 16 CG-nodes, the d-CG scheme 
still faithfully reproduces the phase transition curves and 
fluctuations properties. Since TV^ is largely reduced com- 
pared to N, a considerable speed-up of CPU time, about 
a factor of 40, is realized for N = 16384. A significantly 
higher gain can be expected as the network gets larger, 
allowing the computational study of network size effect 
very affordable. 

Using the d-CG approach with a fixed size of CG- 
networks N'^ = 64, we calculate the dependence of phase- 



transition critical Tc on the network size N, as reported 
in Figl2l^c). It had been well established that the Ising 
model with ferromagnetic interactions on BA scale-free 
network undergoes a phase transition from fcrromag- 
netism to paramagnetism at a critical temperature Tc 
that increases as the logarithm of network size [23|, l29l - 
|3l|. Our result is consistent with the theoretical expres- 
sion dl] that Tc = -^ In N, where (k) = 20 for the 
present case of study. Note that the critical phenomenon 
disappears when N -^ oo. 

For the SIS model there is an epidemic threshold Ac on 
a finite size BA scale- free network network [2J, [S^] ■ Our 
numerical simulation starts from a random configuration 
with about half nodes being infected. After an initial 
transient regime, the system will evolve into a steady 
state with a constant average density of infected nodes. 
The steady density of infected nodes p is computed by 
averaging over at least 50 different initial configurations 
and at least 10 different network realizations with the 
same parameter A. The epidemic threshold Ac occurs at 
p — (absorbing state) if A < Ac and p > (active state) 
if A > Ac [32]. Due to finite size effects, the fluctuation 
can drive the system to the absorbing state, especially in 
the vicinity of Ac. Once the absorbing state is arrived, the 
system will never leave it. Based on the consideration, we 
use, in practice, a nonzero tolerance in p (with the order 
of N~^) as the boundary of the phase transition point. 
Reported in FigHfa) and (b) are the calculated results of 
p and its fluctuation Ap^ = <Jp^) — (p) , respectively, as a 
function of A, obtained by the CG model and the micro- 
scopic model. Again, while the RM scheme fails when the 
relative infection rate parameter A < 0.1, the agreement 
between the d-CG and the microscopic results remains 
excellent. Fixing N'^ ~ 64 in the d-GG approach, the 
resulting Ac shown in FigHKc) is proportional to 1/lnA^, 
also consistent with the theoretical prediction i24j . 

Figure [5] show that the simulation results of the micro- 
scopic and CG levels on scale-free networks with other 
two scaling exponents, 7 = 2 and 7 = 3.5. For larger 
7, Tc becomes lower, while Ac gets larger. It is clearly 
observed that our d-GG approach is still applicable. In 
addition, many other types of networks such as random 
network and small- world network are also tested, and all 
results show that the validity of our d-CG approach does 
not depend on network topology. 



C. Degree correlated netAvorks 

It is worthy noting that the above numerical demon- 
strations are carried out on degree uncorrelated net- 
works. We will show that our d-CG approach is valid 
to reproduce critical behaviors on degree-degree cor- 
related networks as well. It has been witnessed that 
many real networks display different degree-mixing pat- 
terns |33| . To measure the degree of the correlation, in 
Ref.[33| Newman introduced a degree- mixing coefficient: 
Tk = {{kikj) - (fcj) {kj))/{{kf) - (hf), where k^ and kj 




FIG. 5; Comparison of microscopic and CG simulations on 
scale-free networks with the scaling exponent 7 = 2.0 and 
7 — 3.5. Left panel: \m\ and x as a function of T for Ising 
model. Right panel: p and and Ap^ as a function of A for SIS 
model. All networks have fixed size oi N — 1024 and mean 
degree of (k) — 20. 



are the remaining degrees at the two ends of a link and 
(•) means the average over all links, r^ ~ indicates 
that there is no degree correlation, while r^ > (< 0) in- 
dicates that a network is assortatively (disassortatively) 
mixed by degree. Previous studies have revealed that 
degree-mixing pattern plays an important role in dynam- 
ical behaviors on networks, such as percolation 33[, epi- 
demic spreading [3J|, synchronization [3a]- To generate 
different degree-mixing networks, we employ a algorithm 
proposed in [361]. At each elementary step, two links in 
a given network with four different nodes are randomly 
selected. To get an assortative network, the links are 
rewired in such a way that one link connects the two 
nodes with the smaller degrees and the other connects 
the two nodes with the larger degrees. Multiple connec- 
tions are forbidden in this process. Repeat this oper- 
ation until an assortative network is generated without 
changing the node degrees of the original network. Simi- 
larly, a disassortative network can be produced with the 
rewiring operation in the mirror method. We start from 
BA scale-free networks with a neutrally degree-mixing 
pattern, and produce some groups of degree-mixing net- 
works by performing the above algorithm. Figure [5] dis- 
plays the results of Ising model and SIS model for three 
different values of r^ . For each r^ , the simulation results 
of CG models agree well with those of microscopic ones. 
In Ising model, Tc shifts to right and x ^-t Tc becomes 
smaller as rk increases. In SIS model, both Ac and the 
fluctuation of p at Ac decrease with rk ■ 
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FIG. 6: Comparison of microscopic and CG simulations on 
correlated networks for different degree-mixing coefficient r^. 
Left panel: \m\ and x as a function of T for Ising model. Right 
panel: p and and Ap^ as a function of A for SIS model. All 
networks have power-law degree distributions with the scaling 
exponent 7 = 3, and fixed size oi N = 1024 and mean degree 
of (fc> = 20. 



IV. CONCLUSIONS AND DISCUSSION 

In summary, we propose an approach for coarse grain- 
ing the phase transition dynamics on complex networks 
described by stochastic models for both equilibrium 
and nonequilibrium systems. The d-CG approaches via 
degree-based merging scheme are feasible since the re- 
liable microscopic information such as phase transition 
behaviors and fluctuations are preserved. We have veri- 
fied that our d-CG approach supports, in consistent with 
the microscopic models, the equilibrium distributions for 
Ising model and the nonequilibrium dynamical flows for 
SIS model. Stochastic description, as exemplified here, 
is ubiquitously important in the study of phase transi- 
tion dynamics and complex networks for a wide range of 
realistic systems. Moreover, this work also suggests the 
development of other promising CG statistical models 
satisfying CSC. 

It is interesting to compare our d-CG approach with 
heterogeneous mean-field theory (HMFT) that success- 
fully predicts Tc and Ac on heterogeneous scale-free net- 
works [3l|, (sj. Based on the ansatz that nodes with 



the same degree share the same dynamical properties, 
HMFT derives a series of coupled mean-field equations 
for degree-dependent quantities. Recently, Langevin ap- 
proach together with the HMFT has been developed in 
Refs.[37l. |38| . which have confirmed that such an ap- 
proach is responsible for the effect of fiuctuations in a 
finite-size network that often play important roles in the 
vicinity of phase transitions. With regard to the present 
study, we develop a reliable CG simulation approach, 
that not only can correctly predict the critical phenom- 
ena and fluctuations information, but also can be applied 
to diverse networks. Especially, it is shown that the valid- 
ity of our d-CG approach on the application of correlated 
networks; however, in this case the HMFT becomes, in 
general, very difficult to deal with. This is because that 
specific formulation of degree-degree correlation are un- 
known for most of correlated networks, such as P{k'\k), 
that is the conditional probability of a node of degree 
k being connected to a node of degree k'. On the other 
hand, the CSC discussed in this work may provide a solid 
understanding of the physical mechanism behind the ba- 
sic assumption of HMFT. Our analysis here may lead 
to the advancement in efficient and consistent CG ap- 
proaches for dynamics on surfaces and soft lattices. 

Note that the present study is limited to the case of 
quenched networks, that is, the connectivity of networks 
is frozen in time. While for the case of annealed net- 
works, i.e. the networks themselves are dynamical ob- 
jects, our CG approach will encounter some difficulties 
in application. In this case, since the network connec- 
tions are frequently reshuffled, the adjacency matrix of 
the resulting CG-network and CG variables should be 
accordingly updated. This will lead to the very ineffl- 
ciency of our CG approach in simulations. However, an 
important advancement in studying critical phenomena 
on annealed networks has been made in Rcfs.[25, 39] by 
means of finite-size scaling theory. Developing coarse- 
grained simulation methods on annealed networks and 
adaptive networks [40| deserve further investigations. 
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